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RESUMEN 

En la Universidad de Fuerzas Armadas ESPE, desde el 2009, el dector 
Roberto Aguiar viene desarrollando el sistema de computacion CEINCI-LAB para 
el analisis estatico y dinamico de estructuras. Hasta el momenta se ha utilizado el 
modelo bilineal para determinar el desempeno del aislador de Triple Pendulo de 
Friccion. La presente investigacion plantea comparar el modelo bilineal con el tri- 
lineal propuesto por McVitty y Constantinou (2015). Para ello se va a desarrollar un 
nuevo aplicativo en el software Matlab del modelo tri-lineal para un analisis 
respuesta en el tiempo con un grado de libertad e incorporar al CEINCI-LAB. 

Palabras Claves: Aislador de Triple Pendulo de Friccion. Modelo bilineal. Modelo 
tri-lineal. Sistema de computacion CEINCI-LAB. 

ABSTRACT 

At Universidad de las Fuerzas Armadas ESPE, Roberto Aguiar has been 
developing the CEINCI-LAB computing system tar static and dynamic structural 
analysis since 2009. Until now the bilinear model has been used to determine the 
performance of Triple Frictional Pendulum bearing. This paper will compare the 
bilinear model with the tri-linear model proposed by McVitty y Constantinou (2015). 
For this comparison, a new application developed with Matlab software of the tri- 
linear model will be used for the response time analysis with one degree of freedom, 
and will incorporate CEINCI-LAB. 

Keywords: Triple Frictional Pendulum bearing. Bilinear model. Tri-linear model. 
System of CEINCI-LAB computer. 
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1. INTRODUCTION 

El Ecuador se encuentra ubicado dentro del cinturon de fuego del Pacifico, 
que es la zona de mas alto riesgo sismico y volcanico en el mundo, como se pudo 
evidenciar el pasado 16 de abril de 2016, con un sismo de magnitud 7.8. Estos 
eventos sismicos y el avance tecnologico en la construccion de edificaciones con 
aislacion sismica han motivado al Ecuador a no quedarse atras. Actualmente en el 
pais existe una gran cantidad de estructuras con aislacion sismica, se detalla a 
continuacion algunos de ellos: 

Puente “Los Caras”: Este puente se encuentra ubicado al node de la 
provincia de Manabi, une a San Vicente con Bahia de Caraquez, con una longitud 
de 1980 metros. Dispone de 160 aisladores de Triple Pendulo de Friccion, FPT 
(Triple Frictional Pendulum) y los ubicaron entre las bases del puente que se 
encuentran sobre el agua y la via por donde circulan los vehiculos. 



Figura 1 Puente "Los Caras" en Bahia de Caraquez 
Fu e nte : http ://www. tcte levisi o n . co m/ 

Puentes de Esmeraldas: La construccion de los puentes sobre el Estuario 
del rio Esmeraldas permiten el acceso al Puerto Maritimo, al Aeropuerto de Tachina 
y el ingreso a la ciudad de Esmeraldas, en la Figura 2 se observa la ubicacion de 
cada uno, tres de los puentes se encuentran construidos con aisladores de base 
FPT. 



Figura 2 Ubicacion de los Puentes. 
Fuente: Vergara Navarrete (2015) 


Los Puentes construidos con aisladores FPT son: el Puente Sur, el Puente 
Norte y el Puente Norte 1. El Puente Principal de longitud de 434 m se encuentra 
construido con apoyos de neopreno. 
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■ Puente Sur: se encuentra ubicado en el sur de la Isla Prado y en uno de los 
brazos que complementa el Estuario del rio Esmeraldas, tiene una longitud 
de 160 metros, un anchode 18,90 metros y 12 aisladoresde Triple Pendulo 
de friccion, tres por cada estribo y tres por pila, con una capacidad de 
desplazamiento de 23 pulgadas (58,42 cm) +/- 0,3 pulgadas (0,762 cm). 



Figura 3 Puente Sur 
Fuente: Guaygua (2015) 


■ Puente Norte: se encuentra ubicado al node de la Isla Prado y en uno de 
los brazos que complementa el Estuario del rio Esmeraldas, tiene una 
longitud de 120 metros, un ancho de 18,90 metros y 9 aisladores de Triple 
Pendulo de Friccion, con las mismas caracteristicas del Puente Sur. 



Figura 4 Puente Norte 
Fuente: Guaygua (2015) 


■ Puente Norte 1 : se encuentra ubicado node de la Isla Prado cerca del sector 
de la marina, tiene una longitud de 108 metros, un ancho de 18,90 metros 
y 9 aisladores de Triple Pendulo de friccion. 



Figura 5 Puente Norte 1 
Fuente: Vergara Navarrete (2015) 





98 


Pablo Leiva, Roberto Aguiar 


Aeropuerto Internacional Mariscal Sucre: es el nuevo aeropuerto de 
Quito, se encuentra ubicado en la parroquia de Tababela. Con las siguientes 
caracterfsticas operacionales, Quiport (2013): una altitud de 2400 metros sobre el 
nivel del mar, una pista de 41 00 metros de largo, la Torre de Control es de 41 metros 
de altura, el Parqueo de aeronaves tiene 45 puestos, el area para el terminal de 
pasajeros es 38000 metros cuadrados y el terminal de carga es 14000 metros 
cuadrados. En el Puente Peatonal estan colocados 16 aisladores sfsmicos de dos 
tipos debido a las solicitaciones en el diseno. 



Figura 6 Aisladores Sfsmico en el puente peatonal del Aeropuerto de Quito. 

Fuente: Quiport 

Edificio Sede de UNASUR: se encuentra ubicado en la Mitad del Mundo, 
con un area de 20.000 m2 de construction. Es una estructura que debe responder 
en el rango elastico, nucleo con muros acoplados de hormigon armado sobre 
cimentacion continua, dos volados con celosias de acero (38 y 50 metros) y pisos 
aislados en 3 dimensiones al interior de los volados, en la Figura 7 se presenta el 
Edificio Sede de UNASUR, las celosias de aceros y el aislador 3D. 

Criterio de diseno del sistema de aislacion en los volados: 

• Reducir los niveles de aceleracion horizontal y vertical en los pisos de 
oficinas para evitar el dano a componentes no-estructurales y equipos y la 
conmocion de los ocupantes en caso de un terremoto. 

• Mitigar las vibraciones generadas por actividades humanas en los pisos de 
oficinas. 

• Reducir la demanda sfsmica en los elementos de la estructura en celosfa 
que soporta a los volados. 



a) 


b) 
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c) 

Figura 7 a) Perspectiva del Edificio Sede de UNASUR; b) Volados y celosias de 

aceros; c) Aislador 3D 
Fuente: Suarez (2013) 

Centro de Investigation y Posgrado de la Universidad de las Fuerzas 
Armadas - ESPE: es el nuevo Centro de Investigaciones y Posgrado de la UFA- 
ESPE, para la construccion de este proyecto se han requerido de ocho bloques 
estructurales que contaran con 81 aisladores de Triple Pendulo de Friccion y se 
implantara en un area de 23 mil 338 metros cuadrados. En la Figura 8 se observa 
una perspectiva en 3 dimensiones del proyecto, la implantation y la etapa 
constructiva de los cinco bloques estructurales con aislacion sismica. 
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i) j) 


Figura 8 Centro de Investigaciones y Posgrado UFA-ESPE; a) Perspectiva 
Noreste; b) Ubicacion de los bloques estructurales; Fotografias de la etapa 
constructiva hasta octubre del 2016 c) Bloque 1 ; d) Aislador del Bloque 1 ; e) 
Colocacion de aisladores en el Bloque 2; f) Colocacion de aisladores en el Bloque 
6; g) Bloque 7; h) Aislador Bloque 7; i) Bloque 8; j) Aislador Bloque 8. 

Sky Building: es un nuevo proyecto inmobiliario ubicado en la salida node 
del aeropuerto internacional Jose Joaquin de Olmedo, en Guayaquil. Esta 
edification disenada con aisladores sfsmicos es de 10 pisos, con 216 oficinas y 3 
locales comerciales. Los aisladores se colocaron sobre las columnas del subsuelo 
y el bloque de ascensor y gradas esta separado del bloque principal estructural para 
colocar los aisladores sobre la cimentacion, observar Figura 9. 
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a) b) 

Figura 9 a) Aisladores se encuentra sobre columna de subsuelo; b) Aisladores del 
bloque de ascensor y gradas se halla sobre los plintos. 


En la Tabla 1 se resumen las estructuras con aislacion sismica hasta el 
presente ano 2016 en el Ecuador. 


Tabla 1 Estructuras con aislacion sismica en el Ecuador, hasta el ano 2016 


No. 

Estructuras 

Bloques Estructurales 

1 

Centro de Investigacion y Posgrado de 
la UFA-ESPE 

5 

2 

Edificio Sede de UNASUR 

1 

3 

Sky Building 

1 

4 

Aeropuerto Internacional Mariscal Sucre 

1 

5 

Puente “Los Caras” 

1 

6 

Puente de Esmeraldas Norte 

1 

7 

Puente de Esmeraldas Nortel 

1 

8 

Puente de Esmeraldas Sur 

1 

9 

Puente sobre el Rio Chiche 

1 


En el diseno de estructuras con aislacion sismica se requiere que el sistema 
de aislacion sea capaz de responder las deformaciones y cargas correspondiente 
al sismo maxi mo considerado, MCE. Con el objetivo, que la superestructura se 
mantenga en el rango elastico durante el sismo. Por lo tanto, la verificacion del 
desempeno para el MCE se debe realizar analitica y experimentalmente. El criterio 
de esta verificacion se basa en exhibir el peor escenario sismico posible. Por tal 
motivo el presente estudio, busca realizar la verificacion analitica del desempeno 
del aislador de Triple Pendulo de Friccion con el modelo en la fase II de McVitty y 
Constantinou (2015) y comparar con el modelo bilineal mediante un aplicativo en el 
software Matlab. 

2. MODELO DE MCVITTY Y CONSTANTINOU (2015) 

El aislador Triple Pendulo de Friccion (Triple FP) tiene un comportamiento 
no lineal y complejo. A medida que el desplazamiento aumenta, hay multiples 
cambios en la rigidez. Tipicamente a bajos niveles de fuerza / desplazamiento el 
sistema es muy rigido, en comparacion a un sismo de diseno (DBE) en el que hay 
una menor rigidez y un aumento de la amortiguacion, y el comportamiento tambien 
difiere a un sismo maximo considerado (MCE) donde hay una mayor rigidez para 
controlar los desplazamientos. Esto da lugar a un sistema adaptativo que puede ser 
optimizado para diferentes niveles de rendimiento. 
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En Fenz y Constantinou (2008, 2008) y Fadi y Constantinou (2010) 
presentaron un modelo de cinco fases, en base a los 12 parametros geometricos y 
cuatro coeficientes de roce que definen el comportamiento fuerza-desplazamiento 
del aislador triple FP. En este modelo de McVitty y Constantinou (2015), se reduce 
el numero de parametros geometricos a 6 y los coeficientes de roce a 2, observar 
Figura 10, a). 



Figura 10 a) Seccion transversal y b) Vista en code del aislador de Triple FP 
Fuente: McVitty y Constantinou (2015) 


Reduccion de parametros: 

■ Ri = f?4 » /?2 = 

■ M2 = M3 < Ml = |U4 

■ di = d.4 y d .2 = di, 

■ hi = /14 y /l 2 = hz 


Rieff — Ri hi 

„,* _ Rieff j 
ctj — — — — di 

1 Ri l 


( 1 ) 

( 2 ) 


Donde R t es el radio de curvatura; hi y di es la altura y el desplazamiento 
presentado en la Figura 10; R ieff es el radio de curvatura efectivo; d\ es la 
capacidad de desplazamiento; Mi es el coeficiente de friccion. El subfndice /, varia 
de 1 a 4 debido que existen 4 superficies. 


Se adopta estos parametros al comportamiento Fuerza-Desplazamiento y 
se reduce a tres regimen. 

Regimen I 


El desplazamiento se da solamente en las placas 2 y 3. 


0 < u < u* 


u * = 2(mi — R-2)R 2 eff 


(3) 

(4) 
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F = 



u + n 2 W 


(5) 


Donde u es el desplazamiento lateral del aislador; u*maximo 
desplazamiento en la fase I; F es la fuerza lateral aplicada; W es el peso que gravita 
sobre el aislador. 
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3 
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if 









-Regime 1 




Displacement 


Figura 11 Regimen I 
Fuente: McVitty y Constantinou (2015) 


Regimen II 


El movimiento en el aislador interior Mega a los topes y empieza a 
desplazarse en las superficies 1 y 4. 

u* <u < u** ( 6 ) 


u** = u* + 2 di 

F -(z 


(7) 

( 8 ) 


Donde u es el desplazamiento lateral del aislador; u* maxi mo 
desplazamiento en la fase I; u“ maximo desplazamiento en la fase II; F es la fuerza 
lateral aplicada; W es el peso que gravita sobre el aislador. 



Figura 12 Regimen II 
Fuente: McVitty y Constantinou (2015) 


Regimen III 

Este regimen se presenta cuando el sismo es extremadamente fuerte y el 
aislador interior se topa con los topes exteriores, en estas condiciones empieza a 
deslizarse el aislador interior en las superficies 2 y 3. 


cap 


U** <U< U ( 


( 9 ) 
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u cap — 2di* + 2d 2 * (10) 

/ W \ N / W \ (11) 

Donde u es el desplazamiento lateral del aislador; u** maximo 
desplazamiento en la fase II; u cap maximo desplazamiento en la fase III; F es la 
fuerza lateral aplicada; W es el peso que gravita sobre el aislador. 



El analisis sismico de edificaciones con aislacion sfsmica requiere un 
modelo matematico del comportamiento fuerza-desplazamiento del aislador. De tal 
manera, para validar el modelo tri-lineal propuesto por McVitty y Constantinou 
(2015), el aislador no alcanza el Regimen III para desplazamientos MCE. En 
conclusion, el diagrama fuerza-desplazamiento puede ser modelado en el Regimen 
II como se muestra en la Figura 14. La fuerza en el desplazamiento cero viene dada 
por: 


jjtW 


Fi ~ Oi - Fz) 


eff 

R leff. 


w 


( 12 ) 



Figura 14 Modelo Tri-lineal en el Regimen II 
Fuente: McVitty y Constantinou (2015) 
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3. DESCRIPCION DEL APLICATIVO EN MATLAB 

El aplicativo realiza un analisis no lineal de un grado de libertad para la 
verification analitica del desempeno del aislador de T riple Pendulo de Friccion. Este 
analisis se plantea con un modelo tri-lineal propuesto por McVitty y Constantinou 
(2015). 


El punto de fluencia del modelo tri-lineal, se halla con las ecuaciones que 
definen la fuerza y desplazamiento en el regimen I; y el punto ultimo se encuentra 
con la ecuacion que determina la fuerza y desplazamiento para el regimen II. Es un 
modelo conservador ya que se tiene en la realidad como punto ultimo las ecuaciones 
del regimen III. 


En este apartado se describe brevemente los aspectos fundamentales del 
analisis no lineal, para el modelo constitutive tri-lineal. 




Figura 15 Identification de las ramas y puntos de quiebre del Modelo tri-lineal 


A la izquierda de la Figura 15 se definen las 7 ramas que considera el 
modelo, en la rama 1 trabaja exclusivamente en el rango elastico con una rigidez 
Ke que es la pendiente de esta recta, cuando se supera la fuerza de fluencia fy se 
ingresa a la rama 2 denominada post fluencia cuya pendiente es la rigidez kp. En la 
rama dos se va controlando la velocidad, cuando esta Mega a cero se alcanza el 
punto de control 2, que se indica a la derecha de la Figura 15; en la rama 3 se queda 
un momento en reposo hasta que la fuerza lateral supera a la fuerza de friccion, 
2 p 2 VF; ingresando a la rama de descarga definida por 4 con una rigidez Ke. La rama 
5 es similar a la rama 2, la rama 6 a la rama 3 y la rama 7 a la rama 4 pero actuando 
el sismo en sentido contrario. En analisis no lineal, es muy importante la 
determination de los puntos de quiebre, en forma exacta y para lograr este objetivo 
se disminuye el incremento de tiempo At. Revisar Aguiar (2012) y Freire y Ortiz 
(2014). 


3.1. Metodo p de newmark 

El marco teorico del metodo /? de Newmark, se halla descrito en Aguiar 
(2012) y Freire y Ortiz (2014). Aqui se presenta el procedimiento computacional 
para el analisis no lineal. Partiendo de que inicialmente se halla en la rama 1 ; en 
este caso y al inicio del programa para t = 0, las condiciones iniciales son nulas. 
Con esta introduction el procedimiento de calculo es el siguiente. 
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i. Definir los datos del aislador y acelerograma: 

Geometria del aislador ; Coeficientes de Friction-, 
m = masa del aislador + masa actuante sobre el aislador ; W = peso total’, 
g = Acelerograma ; p = g *m 

ii. Se determina el caso especial que se va a utilizar. 

-> Metodo de la aceleracion constante, y = 1/2 y y? = 1/4 

-> Metodo de la aceleracion lineal, y = 1/2 y /? = 1/6 

-> Metodo de la aceleracion escalonada, y = 1/2 y /? = 1/8 

iii. Se define las Condiciones iniciales: 

q 0 = 0 m q 0 = 0 m/s g 0 = 0 tti/s 2 
Donde q t es el desplazamiento; eg es la velocidad; eg es la aceleracion. 
Calculos que se realizan para cada iteracion i: 

iv. Se determina el intervalo de tiempo At : 

A t = t i+1 -t ; (13) 

v. Se calcula el amortiguamiento: 

Los valores de k t depende de la rama en que se encuentra, observar 
Figura 15 y f se detalla en el literal 3.2. 

c = 2 (14) 

Donde /q es la rigidez elastica o plastica; f es el factor de 
amortiguamiento; c es el amortiguamiento para un sistema de un grado de 
libertad. 

vi. Se obtiene la rigidez equivalente y el intervalo de carga: 


1 , y 

a = P* At* m+ P* C 

(15) 

fo = 2*/?* m + At * G/? O*' 

(16) 

y 1 

k; = k:+ n *c + n A ,*m 

' ' P*A t [3 * At 2 

(17) 

A Pi = Api+ a cg + b q. 

(18) 


Donde fc ; es la rigidez equivalente; A p t es el intervalo de carga. 


vii. Se calcula el intervalo de desplazamiento y velocidad: 
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A 

4,1 = T 

(19) 


(20) 

Finalmente se obtiene la respuesta de desplazamiento, velocidad, la fuerza 
lateral y aceleracion del aislador en cada instante de tiempo: 

<7i+ 1 = + <7; 

(21) 

q l+ i = A q + cii 

(22) 


ix. La aceleracion se calcula preferible mediante la ecuacion diferencial de 
movimiento: 

.. _ Pi+i ~ c * <7i+i ~ (fs) i+1 (23) 


Se destaca que: 

(fs ) i+ 1 = (fs) t + ki A q { (24) 

Donde fs es la fuerza lateral asociada a la rigidez, que va cambiando de 
acuerdo a la rama en que se encuentra el modelo. 

3.2. Factor de amortiguamiento (£) 

En la fase I el aislador empieza con un factor de amortiguamiento igual a 
0.05, debido que el aislador empieza trabajando con una rigidez alta y un 
amortiguamiento bajo. Cuando pasa a la fase II, se calcula en cada instante el factor 
de amortiguamiento con la siguiente ecuacion: 

E = (25) 

? 4tt£ s 



Figura 16 Calculo de la Energfa Disipada y Energfa Elastica. 
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En la Figura 16 el eje de las abscisas y ordenadas divide en 4 areas iguales 
(F * £>£>) , y para encontrar la Energia Disipada por cada ciclo se resta las areas 
externas, Ai , A2 y A3; como resultado la siguiente ecuacion, se describe en Gutierrez 
(2016): 

E d = 4 (F * £>£>) - 2A t - 2A 2 - 2A 3 (26) 

Simplificando los terminos: 

E d = 4 * (q ± * \l 2 * W — q ± * F + DD * 2 * ^ * W ) (27) 

Metodo de Calculo 

i. Se determina en cada instante de tiempo el desplazamiento, DD 

ii. Con el desplazamiento se calcula su Fuerza, F. 

F = K 2 *DD + h*W (28) 

iii. Se calcula la Energia Disipada 

E D =4*{q 1 * [ i2*W-q 1 *F + DD * 2 * \a 1 *W') (27) 

iv. La Energia Elastica 

F * DD (29) 


v. Zeda para cada instante. 

, = _^ d _ (25) 

? 4t zE s 


4. ENSAYO DEL AISLADOR FPT8833/1 2-1 2/8-5 

El aislador FPT8833/12-12/8-5 del nuevo Centro de Investigacion y Posgrado 
de la ESPE, fue ensayado bajo una funcion armonica que se presenta en la Figura 
17. El periodo promedio de los bloques con aisladores sfsmicos para el DBE es Td 
= 3s. Por lo tanto, para el ensayo del aislador se utilizo un frecuencia io n =2Tr/TD. 
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Figura 17 Lafuncion armonica a) Desplazamiento-Tiempo; b) Velocidad-Tiempo; 

c) Aceleracion-Tiempo 


Se observa en la Figura 18 la constitutiva del ensayo experimental que es 
ajustada a un modelo matematico para determinar los coeficientes de friccion e 
imagenes del ensayo. En la Figura 19 se presenta la constitutiva analitica del 
modelo tri-lineal propuesto por McVitty y Constantinou (2015) aplicado la misma 
carga del ensayo experimental en el programa de Matlab. 




c) 

Figura 18 a) Instalacion del aislador en el ensayo experimental; b) Aplicacion de 
la carga al aislador c) Constitutiva ajustada del ensayo experimental bajo tres 

ciclos de carga. 

Fuente: EPS 
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Figura 19 Constitutiva obtenida en este estudio con el modelo tri-lineal de McVitty 

y Constantinou (2015) 

En el ensayo experimental empieza con una fuerza en el desplazamiento 
igual a cero, para obtener en los 3 ciclos desplazamientos semejantes. En cambio, 
en el aplicativo de Matlab se observa que en el ciclo 1 se obtiene un desplazamiento 
menor porque parte de la fase 1 y los dos ciclos siguientes dan el desplazamiento 
final. 


El desplazamiento promedio del ensayo experimental en los 3 ciclos es de 
342.3 mm. En cambio, mediante el aplicativo en Matlab se obtuvo un 
desplazamiento de 334.3 mm. Este dato permite verificar la programacion del 
aplicativo del modelo matematico que plantea McVitty y Constantinou (2015) y se lo 
puede utilizar para realizar un analisis no linear respuesta en el tiempo del aislador 
de Triple Pendulo de Friccion en un sistema de 1 grado de libertad. 

4.1. Coeficientes de Friccion 

Para determinar los coeficientes de friccion del ensayo experimental se 
ajusta las graficas a un modelo matematico, en este estudio se va explicar mediante 
el modelo de la fase II de McVitty y Constantinou (2015). A continuation se 
determina el proceso de calculo: 



Figura 20 Modelo ajustado del ensayo experimental para determinar los 
coeficientes de friccion 

Del modelo ajustado se extrae los dos valores que se muestra en la Figura 
20 a la izquierda y se iguala con el modelo de la fase II de McVitty y Constantinou 
(2015). Primero se obtiene g2 (5.36=2g2W) y para pi se despeja del valor de la 
fuerza en el desplazamiento cero (14.56=pW), se conoce que: 
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HW = 


Ai - (A i - Az) 


file//. 


W 


( 12 ) 


5. PELIGROSIDAD SISMICA Y ESPECTROS DE ANALISIS 

Quito se halla sobre cuatro fallas geologicas inversas, denominadas: 
Puengasi, llumbisi-La Bota, Carcelen el Inca y Bellavista-Catequilla, cuyas trazas y 
pianos de falla se presentan en la Figura 21 . 



Figura 21 Modelos de Fallas ciegas de Quito, trazas y pianos de falla. 
Fuente: Rivas et al. (2014) 

En el estudio de “Espectros de control para la ciudad de Quito-Ecuador” 
las mayores ordenadas espectrales para las zonas: 

• Sur, centra sur y centra de Quito: se hallaron ante un sismo de magnitud 6.4 
registrado en el centra del piano de falla del segmento de Puengasi. 

• Centro norte y node de Quito: se hallaron ante un sismo de magnitud 6.2, en la 
falla de llumbisi-La Bota. 

La falla de Puengasi tiene un periodo de recurrencia de 1224 y 2190 anos 
para que se produzca un sismo Mw=6.4, valores hallados con modelo de Gutenberg 
y Richter Modificado. Rivas et al. (2014). 



Figura 22 Espectros encontrados en perfil tipo C, para el sector Sur de Quito. 

Fuente: Aguiar (2015) 


En la Figura 22 se presenta los espectros de analisis para la falla de 
Puengasi medidos en un perfil de suelo tipo “C” en el Sur de Quito, para niveles de 
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confianza del 50% (Ifnea continua), 84% (Ifnea discontinua) y el espectro hallado 
con la NEC-15, el DBE. 

En el perfodo de 1 segundo, los espectros hallados con un nivel de 
confianza del 84% superan la aceleracion espectral encontrada con el NEC-15. Por 
tal motivo, para encontrar el sismo MCE a partir del sismo DBE se busca un factor, 
que es la division hallada con el modelo de Campbell y Borzognia (2014) (Ifnea 
naranja entrecortada) para la aceleracion del NEC-15 y su valor es de 1.4. 

5.1. Acelerogramas seleccionados 

Los registros que se utilizaran deben asemejarse a las caracterfsticas 
geologicas de Quito, con las siguientes consideraciones que proponen Aguiar y 
Morales (2016): 

• Los registros son de fallas cercanas a las fallas. 

• Los registros del sismo de Gazly de 1976, provienen de una region tipo 
“backward-directivity”. Los restantes registros son de regiones “forward- 
directivity”. 

• La magnitud de los registros se encuentra entre 6.7 y 7.1 . 

• Las distancias a la fuente Rrup se encuentran entre 3 y 1 2 km, 

• Los registros fueron obtenidos en suelos tipo “C” y “D” 

• Los registros de Gazli de 1976 y el de Loma Prieta de 1989, son considerados 
como registros estandar con alto contenido de frecuencia. 

• Los sismos de Northridge de 1 994 y Kobe de 1 995 son sismos impulsivos y son 
considerados por su cercanfa a la falla. 

• Se incluye el sismo de Turqufa de 1999 porque fue altamente destructive. 


Tabla 2 

Registro de sismos seleccionados 

Codigo 

Nombre del Sismo 

Estacion 

M w 

R(Km) 

Tipo de 
Suelo 

EQ1 

1976 Gazli, USSR 

Karakyr 

6.80 

5.46 

C 

EQ2 

1989 Loma Prieta 

LGPC 

6.93 

3.88 

C 

EQ3 

1989 Loma Prieta 

Saratoga-W 
Valley Coll 

6.93 

9.31 

c 

EQ4 

1994 Northridge 

Jensen Filter 
Plan 

6.69 

5.43 

c 

EQ5 

1994 Northridge 

Sylmar-Coverter 
Sta. East 

6.69 

5.19 

c 

EQ6 

1995 Kobe, Japon 

Takarazuka 

6.90 

3.00 

D 

EQ7 

1999 Duzce, 
Turquia 

Bolu 

7.14 

12.41 

D 


5.2. Espectros de respuesta elasticos 

El espectro de diseno denominado DBE, que corresponde al NEC-15 y el 
maximo considerado MCE (observar la Figura 23) que se obtiene al multiplicar el 
DBE por 1 .4, dicho factor se determino en el apartado anterior. 
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En la Figura 24 estan los espectros obtenidos de cada acelerograma en la 
componente normal y paralela a la falla para un factor de amortiguamiento del 5%; 
la linea roja es el promedio de todos los espectros, 50% de nivel de confianza. 


Como se puede observar en la Figura 25 el espectro MCE esta por encima 
del espectro promedio de los registros seleccionados, por lo tanto se debe escalar 
al MCE los registros para los periodos de interes. 



Figura 23 Espectros elasticos de Diseno DBE y maximos considerados MCE 
para Quito en un suelo tipo "C" 



Periodo (s) 


Figura 24 Espectros elasticos de los registros considerados para el Sur de Quito 



Figura 25 Espectros de diseno DBE, maxi mo considerado MCE y promedio de los 

registros seleccionados 

5.3. Factores de escalamiento 

Para el escalamiento de sismos el rango de los periodos es 0.5TDa 1 .25 Tm. El 
periodo de la estructura para el DBE es Td y para el MCE es Tm. Para los cincos 
bloques estructurales de la UFA-ESPE el Td =3s y el Tm =3. 4s. De tal manera el 
rango de escalonamiento es de 1 .5 y 4.25. En la Tabla 3 se presenta el factor de 
escalonamiento para cada registro. 
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Tabla 3 Factor de Escalonamiento mediante el metodo del Promedio de Pesos 

Ponderados. 


FACTOR DE ESCALONAMIENTO 


EQ1 

EQ2 

EQ3 

EQ4 

EQ5 

EQ6 

EQ7 

Sa MCE (t=1) 

0.941 

0.941 

0.941 

0.941 

0.941 

0.941 

0.941 

Sa FN (t=1) 

0.815 

0.873 

0.710 

0.512 

0.732 

0.933 

0.729 

Sa MCE (t=2) 

0.471 

0.471 

0.471 

0.471 

0.471 

0.471 

0.471 

Sa FN (t=2) 

0.198 

0.561 

0.314 

0.411 

0.429 

0.396 

0.259 

Sa MCE (t=3) 

0.314 

0.314 

0.314 

0.314 

0.314 

0.314 

0.314 

Sa FN (t=3) 

0.161 

0.436 

0.176 

0.421 

0.394 

0.148 

0.146 

Sa MCE (t=4) 

0.235 

0.235 

0.235 

0.235 

0.235 

0.235 

0.235 

Sa FN (t=4) 

0.191 

0.327 

0.078 

0.195 

0.186 

0.072 

0.084 

Fj 

1.375 

0.869 

1.461 

1.128 

1.084 

1.144 

1.527 

FI 

1.200 

1.200 

1.200 

1.200 

1.200 

1.200 

1.200 

FI Fj 

1.650 

1.043 

1.754 

1.353 

1.301 

1.373 

1.833 


En la Figura 26 se representa con Ifnea roja el promedio de los registros 
multiplicados por el factor de escalonamiento respectivo, como resultado se tiene 
que las ordenadas esten sobre el MCE. 



0 0.5 1 1.5 2 2.5 3 3.5 4 

Periodo (s) 

Figura 26 Espectro de Diseno DBE, maximo considerado MCE y promedio de los 

7 sismos FN 

6. DESEMPENO DEL AISLADOR FPT8833/1 2-1 2/8-5 


El aislador FPT8833/12-12/8-5 del nuevo Centro de Investigation y 
Posgrado de la ESPE se va a determinar su desempeno con el modelo tri-lineal de 
McVitty y Constantinou (2015) en el sector Sur de Quito para un suelo tipo “C” para 
comparar con los resultados obtenidos en Aguiar y Morales (2016). Esta 
comparacion permite evaluar el modelo tri-lineal con respecto al bilineal, en la Figura 
27 se presenta ambos modelos. 

La geometna y coeficientes de friction del aislador estan descritos en 
Constantinou et al. (2016) 
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Figura 27 Constitutivas utilizadas a) Modelo Tri-lineal; b) Modelo Bilineal 

6.1. Resultados 


Para determinar el desempeno del aislador 8833/12-12/8-5, se realizo 
mediante el modelo de tres fases de McVitty y Constantinou (2015). Los coeficientes 
de friccion se consideraron sus valores nominales y una carga de 268 Toneladas. 

A los acelerogramas seleccionados se les multiplico por el factor de 
escalamiento respectivamente para encontrar la respuesta en el tiempo del aislador 
mediante el modelo tri-lineal (color magenta) y bilineal (color azul), estos resultados 
se presentan en la Figura 28. 

Cabe indicar que el valor del ^ se calcula en cada instante para simular el 
comportamiento no lineal y no trabajar unicamente con el ^efectivo, como se lo 
realizo en el modelo bilineal. 


Adicionalmente se presenta una grafica comparativa del desplazamiento en 
el tiempo de cada registro mediante los dos modelos, para comparar sus 
desplazamientos, observar en la Figura 29. 


Gazli 1 976 Karakyr FN 


Loma Prieta 1989 LGPC FN 



Loma Prieta 1989 Saratoga FN Northridge 1994 Jensen FN 
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Northridge 1 994 Sylmar FN Kobe 1 995 Takarazuka FN 



Turquia 1999 Bolu FN 



Figura 28 Respuesta en el Tiempo del aislador con los 7 registros FN mediante el 
comportamiento (A) Tri-lineal y (B) Bilineal. 


Gazli 1976 Karakyr FN Loma Prieta 1989 LGPC FN 



Loma Prieta 1989 Saratoga FN 


Northridge 1994 Sylmar FN 


Northridge 1994 Jensen FN 


Kobe 1995 Takarazuka FN 



Turquia 1999 Bolu FN 



Figura 29 Grafica comparativa del desplazamiento en el tiempo mediante el 
modelo tri-lineal y bilineal. 
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Mediante el nuevo aplicativo en Matlab se observa mejores resultados 
porque el programa esta ejecutandose mientras se tiene las mayores aceleraciones 
del registro y esto permite verificar analiticamente el desplazamiento maximo 
esperado. Pero se aprecia que deja de correr el programa cuando las aceleraciones 
tienden a cero debido que el aislador es un sistema adaptativo y para bajas 
condiciones de fuerzas es muy rigido. 

En la Figura 29 se aprecia claramente que el desplazamiento en el tiempo 
tiene el mismo comportamiento en ambos modelos, que era el resultado que se 
esperaba para corroborar los modelos. 

En la Tabla 4 esta detallado los desplazamiento maximos esperados en 
cada sismo mediante los dos modelos y el mayor desplazamiento obtenido con el 
modelos tri-lineal es de 46 cm y con el bilineal es 46.4 cm. El sismo de Loma Prieta 
es el que genera mayores desplazamientos. La capacidad de desplazamiento 
maximo del aislador en la fase II es 53.3 cm, de tal manera se verifica analiticamente 
que el aislador FPT8833/1 2-1 2/8-5 trabaja perfectamente bajo las condiciones 
geologicas del sector Sur de Quito en un suelo tipo “C”. 


Tabla 4 Validation de los desplazamientos en la Fase II mediante el modelo Tri- 
lineal y Bilineal en los 7 registros. 


Codigo 

Nombre del 
Sismo 

Estacion 

Desplazamiento 

Max(mm) 

Tri-lineal Bilineal 

Diferencia 

desplazamientos 

(%) 

EQ1 

1976 Gazli, 
USSR 

Karakyr 

258.9874 

258.5895 

-0.15 

EQ2 

1989 Loma 
Prieta 

LGPC 

460.7535 

464.3997 

0.79 

EQ3 

1989 Loma 
Prieta 

Saratoga-W 
Valley Coll 

317.7102 

312.78 

-1.55 

EQ4 

1994 

Northridge 

Jensen 
Filter Plan 

426.3481 

445.1043 

4.40 

EQ5 

1994 

Northridge 

Sylmar- 
Coverter 
Sta. East 

415.9048 

431.5528 

3.76 

EQ6 

1995 Kobe, 
Japon 

Takarazuka 

134.534 

138.6055 

3.03 

EQ7 

1999 Duzce, 
Turquia 

Bolu 

258.1822 

239.668 

-7.17 


7. CONCLUSIONES 

El comportamiento del aislador de Triple Pendulo de Friccion es no lineal y 
complejo. Con el modelo matematico que propone McVitty y Constantinou (2015) 
permite verificar analiticamente el desempeno del aislador, como se pudo evidenciar 
en el ensayo experimental se obtuvo un desplazamiento de 34.2 cm, y 
analiticamente un desplazamiento de 33.4 cm equivalente a un error del 3%. Con 
este resultado se verifica la utilization para un analisis sismico el modelo tri -lineal 
de la fase II de McVitty y Constantinou (2015). 
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Comparar los resultados del modelo tri-lineal y bilineal para el analisis no 
lineal de un sistema de un grado de libertad permite afirmar que ambos modelos 
son utiles para determinar el desempeno del aislador. En la Tabla 4 se presenta los 
desplazamientos maximos para cada registro y la diferencia en desplazamiento 
entre ambos modelo, la mas alta se obtuvo en el sismo de Duzce 1999 del 7%, esta 
cifra mas el 3% de error obtenida en el ensayo experimental es 10%. Por lo tanto, 
para un prediseno si se recomienda trabajar con el modelo bilineal pero para el 
diseno final se deberfa trabajar con el modelo tri-lineal de la fase II. 

Al aislador FPT8833/1 2-1 2/8-5 utilizado en el nuevo Centro de Investigation 
y Posgrado de la UFA-ESPE se determino su desempeno bajo la demanda sfsmica 
encontrada en el presente estudio para el sector Sur de Quito en un suelo tipo “C”, 
con un desplazamiento maximo de 46 cm en el sismo de Loma Prieta 1989, 
previamente escalado al espectro maximo considerado. Se concluye que el uso del 
aislador FPT8833/1 2-1 2/8-5 se comporta favorablemente porque la capacidad de 
desplazamiento es mayor a la demanda sfsmica calculada. 
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